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ABSTRACT 

We present the most precise to date orbital and physical parameters of the well known 
short period (P = 5.975 d), eccentric (e = 0.3) double-lined spectroscopic binary BY 
Draconis, a prototype of a class of late-type, active, spotted flare stars. We calculate 

■ the full spectroscopic/astrometric orbital solution by combining our precise radial ve- 
Q\ ' locities (RVs) and the archival astrometric measurements from the Palomar Testbed 
t/"") i Interferometer (PTI). The RVs were derived based on the high resolution echelle spec- 
tra taken between 2004 and 2008 with the Keck I/HIRES, Shane/CAT/HamSpec 

*0 | and TNG/SARG telescopes/spectrographs using our novel iodine-cell technique for 

. double-lined binary stars. The RVs and available PTI astrometric data spanning over 

8 years allow us to reach 0.2-0.5% level of precision in Msin 3 i and the parallax but 
the geometry of the orbit (i ~ 154°) hampers the absolute mass precision to 3.3%, 
which is still an order of magnitude better than for previous studies. We compare 
our results with a set of Yonsei-Yale theoretical stellar isochrones and conclude that 
BY Dra is probably a main sequence system more metal-rich than the Sun. Using 
>^ ■ the orbital inclination and the available rotational velocities of the components, we 

also conclude that the rotational axes of the components are likely misaligned with 
the orbital angular momentum. Given BY Dra's main sequence status, late spectral 
type and the relatively short orbital period, its high orbital eccentricity and probable 
spin-orbit misalignment are not in agreement with the tidal theory. This disagreement 
may possibly be explained by smaller rotational velocities of the components and the 
presence of a substellar mass companion to BY Dra AB. 

Key words: binaries: spectroscopic - binaries: visual - stars: fundamental parameters 
- stars: individual (BY Dra) - techniques: interferometric - techniques: radial velocities 



1 INTRODUCTION 



Regular studies of BY Dra (Gl 719, HD 234677, HIP 91009, 
NLTT 46684, BD+51 2402) started in mid 40's when Munch 
noted the calcium H and K lines to be in emiss i on (Miinchl 
1 19441 ). This fact was later confirmed bv lPopperl (1 19531 ) who 
also noted strong emission Balmer series lines. In one of 
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his spectrograms the emission was particularly strong which 
led to a conclusion th at BY Dra ma y be a member of a new 
group of flare stars (|Poppeij fl953l). Photometric monitor- 
ing was then carried out (e.g. Masani et al.|[l955l ) but the 
variability was not confirmed until 1966 when Chugainov 
obtained a quasi-sinusoidal light curve with an amplitud e 
of 0.23 mag and a period of 3.826 d (I Chugainovl 1 1 9661) ■ 
later interpreted as a rotation of a spotted star 1 Krzemiriskil 
1969). No flares were then observed. The first photomet- 
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ric flares were reported by ICristaldi fc Rodonol (|l968l ) who 
later observed twelve flares that occurred between July 1967 
and Ju ly 1970 (|Cristaldi fc Rodond[l970l . ll97ll ). lKrzemihskil 
(1969) confirmed the sinusoidal variability with a ~3.826 d 
period and noted the variation of its amplitude. Many sub- 
sequent studies of BY Dra's variability have been carried 
out and the most up-to-date value of the rotational period 
P ro t = 3.8285 d is given by IPettersen. Olah fe Sandmannl 
l| 19921 ) as an average period for their entire 1965-1989 data 
set. 

The double-line d spect roscopic nature was revealed by 
iKrzemihski fc KraftJ (|l967l ). They announced a period of 
5.981 d but neve r published the i r full orbital solution. It 
was done later by iBopp fc Evansl (|l973h on the basis of 23 
spectra taken between June 1966 and July 1971 at the Hale 
and McDonald obser vatories, 15 of w h ich sh owed unblended 
Ca II H and K lines. IBopp fc Evans! (| 19731) also performed 
an analysis of the spots on the surface of BY Dra and esti- 
mated the rotational velocity of the primary (spotted) com- 
ponent to be ~5 km s _1 and the (rotational) inclination 
to be ~30 deg. Since then BY Dra became a prototype of 
a new class of stars characterized by a late type, bright- 
ness variation caused by spots, rapid rotation and strong 
emission in H and K lines. The short orbital period also 
seems to be a characteristic for most of BY Dra-type stars 
l|Bopp. Noah fc Klimkelll980l ). 

The or bi tal s olution was later improved by 
IVogt fc Fekel (Il979h on th e basi s of high-resolution 



reticon spectra. Vogt fc Fekell l| 19791) also found the pro- 
jected rotational velocity of the primary to be 8.5 km s~ 
under the assumption of the rotational inclination being 
the same as the orbital one (spin-orbit alignment). They 
estimated the radius of the primary to be greater than 0.9 
R,0 which led to a conclusion that BY Dra was a pre-main- 
sequence system. This conclusion was supported by the 
large brightness ratio despite the mass ratio being close to 
1 (q — 0.98), the Barnes-Evans visual surface brightness re- 
lation l|Barnes. Evans fc Moffetlll978l ) and the inequality of 
the rotational and orbital periods. However, the assumption 
of the spin-orbit alignment in close binary systems was later 
critici zed in several works, e.g. Gl ebocki fc Stawikowskil 
(1 19951 ). The orbital parameters as well as the value of the 
projected rotational velocit ies for both componen ts were 
shortly after improved by iLucke fc Mayor] (|l980l ). They 
used new measurements from CORAVEL and obtained 
rotational velocities wisini = 8.05 and U2sini = 7.42 
km s _1 and the mass ratio q = 0.89 sign i ficantl y more 
different from 1 than that of I Vogt fc Fekell (|l979l ). They 
also estimated the magnitude difference (1.15 mag) and 
the primary's radius (1.2 - 1.4 R©) but noted that a 
higher macroturbulence velocity would reduce the radius 
estimation by a factor of 2. 

The most up-to-date spectroscopic orbital solution 
was given together w ith the first astrometric solution by 
iBoden fc Lane! (2001). They combined the archival RV mea- 
surements with the visibility based (V 2 ) astrometric mea- 
surements obta2ned_w^the_Palomar Testbed Interferome- 
ter (PTI; IColavita et al.|[l999l ) in 1999. Their orbital incli- 
nation (152 deg, retrograde orbit) agr ees with the first esti - 
mations of the rotational inclination (IBopp fc Evans)ll973l ) 
but not with the later ones (|GIebocki fc Stawikowskill 19951) . 



cal multiple system. A common prop er-motion companion 
was found bv lZuckerman etakl |l997l ) about 16.7 arc-sec to 
the northeast of the primary. From the visual and infrared 
photometry of BY Dra C they also deduced that this com- 
ponent is a normal M5 dwarf at least 3 x 10 s yr old which 
makes the pre-main-sequence nature of BY Dra less proba- 
ble. Yet another putative component is repor ted in the H ip- 
parcos Double and Multiple System Annex (|ESAlll997l ). A 
photocentric circular orbital solution with a period o f 114 d 
and 113 deg inclination is reported. IBoden fc Land (|200ll ) 
however demonstrated that this is an improbable solution 
since the 4-th body would produce significant perturbations 
to the BY Dra AB radial velocities but no such periodicity 
is seen in the archival RVs. 

We spectroscopically observed BY Dra over the years 
2004-2008 using a combination of high resolution echelle 
spectrographs HIRES (10-m Keck I), SARG (3.5-m TNG) 
and HamSpec (3-m Shane telescope) a s a part of our ongo- 
ing R V search for circumbinary planets (|Konacki et al.ll2009l . 
2010). Even though we knew that BY Dra was too variable 
to allow us reach an RV precision sufficient to detect plan- 
ets, it was nevertheless observed to make use of an extensive 
and publicly available set of PTI V 2 measurements spanning 
now over 8 years. 

In this paper we present a new orbital solution and the 
orbital and physical parameters of the BY Dra AB binary, 
derived wit h a precision of o ver a n order of magnitude bet- 
ter than bv lBoden fc Lanel (|200lT ) . Thanks to our superior 
iodine-cell-based radial velocities and the full set of PTI vis- 
ibilities we are able to put strong constraints on the nature 
of the system. In Sections 2 and 3 we present the data — 
V 2 s and RVs. In Section 4 we describe their modeling. The 
results of our data modeling are presented in Section 5 and 
the state of BY Dra is then discussed in Section 6. 



2 VISIBILITIES 

Often the main observable in the interferometric observa- 
tions at optical or infrared wavelength is the normalized am- 
plitude of the coherence function - a fringe pattern contrast, 
commonly known as the visibility (squared, V 2 ) of th e inter- 
ferom etric fringes, calculated by definition as follows |Bodenl 
Il999h : 



y2 / J-max Imin \ 

V Imax H~ f-min J 



(i) 



where I m ax and Imin are the maximum and minimum in- 
tensity of the fringe pattern respectively. For a given object 
the observed V 2 depends on its morphology and the pro- 
jected baseline vector of a two-aperture interferometer B± 
onto a plane tangent to the sky. For binaries, V 2 varies also 
due to the orbital motion of the components. In the case of 
a binary, approximated by two uniform dis ks, the squar ed 
visibility can be modeled as follows (see e.g. lBodenlll999l ): 



binary 



V? + r 2 Vj + grViVg cqs(27tB_l ■ As/A) 
(1 + r) 2 



(2) 



Finally it is worth noting that BY Dra is a hierarchi- 



where Vi,a are the visibilities of uniform disks (components) 
of the angular diameters 6*i and 8i and are calculated as 
follows: 
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Table 1. Absolute values of the radial velocities of BY Dra with their errors and the best-fit O — Cs. The formal error is denoted with 
cr and the adopted final error with e. The subscript "1" is for the primary and "2" for the secondary. K/H denotes the measurements 
from Keck 1/HIRES, T/S from TNG/SARG and S/H from Shane/HamSpec. 
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:i = 1,2;B, = IIS, 



(3) 



where r is the brightness ratio at the observing wavelength A, 
Ji(x) is the first order Bessel function and As = (Aq, AS) 
is the separation vector between the primary and the sec- 
ondary in the plane tangent to the sky. This vector is related 
to the Keplerian orbital elements, orbital period P, eccen- 
tric anomaly E (from the Kepler equ ation E — e sin E = M ) , 
and the parallax k in the usual way (|van de Kamdfl967h . 

A visibility measurement needs to be calibrated by ob- 
serving at least one calibration source before or after a target 
observation. The calibrator is typically a single star with a 
known diameter and its visibility V c 2 ai is given by Relation 
[31 The correction factor / which should be applied to the 
observed target V 2 is simply the ratio / = V c 2 a ,/V c 2 ai _ meas 
where V 2 al _ meas is the measured calibrator visibility. The 
"true" target visibility is then 

Vtrue J '^measured' (4) 

Uncalibrated visibilities of BY Dra were extracted from 
the NASA Exoplanet Science Institute (NExSci) database 
of the PTI measurement These measurements were made 
in K (2.2 fim) and H (1.6 (im) bands. They were calibrated 
using the standard tools provided by NExSci (getCal and 
wbCalib). As the calibration objects we used HD 177196 
(A7V, V = 5.0 mag, K = 4.5 mag, diameter 9 = 0.42 mas, 
6.6 deg from BY Dra) and HD 185395 (F4V, V = 4.5 mag, 



1 https://nexsciweb.ipac.caltcch.edu/pti- 
archive/secure/main.jsp 



K = 3.5 mag, 9 = 0.73 mas, 9.9 deg) as in Bo den fc Land 
|200ll ). We do not list these measurements as they can be 
easily obtained using the NExSci database and tools. 



3 RADIAL VELOCITIES 

Our high-resolution echelle spectra of BY Dra were obtained 
during 17 nights between September 2004 and November 
2008. We collected 24 spectra using Keck I/HIRES (K/H, 15 
spectra), TNG/SARG (T/S, 4) and Shane/HamSpec (S/H, 
5) telescopes/spectrographs. Our spectra have the resolu- 
tions R ~ 67 000 for K/H, 86 000 for T/S and 60 000 for 
S/H. The typical signal to noise ratio (SNR) per collapsed 
pixel at 550 nm was ~250 for K/H, ~90 for T/S and ~60 
for S/H. The basic reduction (bias, dark, flatfield, scattered 
light subtraction) was done with the CCDRED and echelle 
packages from IRAIQ The wavelength solution and radial 
velocities were obtained with our nove l procedure based on 
the io dine cell technique (|Kon acki 2009; Ko nacki et al.l 120091 . 
2010). This procedure employs a tomographic disentangling 
of the component spectra of double-lined spectroscopic bi- 
naries (SB2s) implemented through a maximum entropy 
method and the two-dimensional cross-correlation technique 



2 I RAF is written and supported by the IRAF programming group 
at the National Optical Astronomy Observatories (NOAO) in 
Tucson, AZ. NOAO is operated by the Association of Uni- 
versities for Research in Astronomy (AURA), Inc. under co- 
operative agreement with the National Science Foundation. 
http://iraf.noao.edu/ 
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TODCOR l|Zucker fc M azch 1991) using syn thetic spectra de- 
rived with ATLAS 9 and ATLAS 12 codes (|Kuruc z 1995) a.s 
templates for the first approximation of the RVs. With this 
approach it is possible to re ach up to 2 m s -1 p recision in 
RVs for components of SB2s |Konacki et al . 2009) but in the 
case of BY Dra the precision is hampered by the activity of 
the star (presence of spots) and the relatively rapid rotation 
of both components. 

In Table \T\ we list our RV measurements together with 
their uncertainties and the best-fit O — Cs. The formal er- 
rors, a, were calculated from the scatter between orders and 
predominantly reflect a high SNR of our spectra. The formal 
errors underestimate the true RV scatter (due to activity) 
and the resulting reduced \ 2 of the spectroscopic orbital fit 
was much larger than 1. Hence to obtain a conservative esti- 
mation of the parameters' errors (and the reduced \ 2 close 
to 1) we added in quadrature a systematic error a 3ys of 150 
m s _1 . Let us note that spots can easily induce RV variations 
at the level of a few hundreds of m s~ so the RV variability 
of BY Dra is not surprising (see e.g. iHelminiak fc Konackil 
l201ll ; IHelminiak et alj l201lh . We also had to adopt small 
shifts between each data set as is explained in lKonacki et al.l 
(2010). The best fit values of the shifts can be found in Ta- 
ble [3J in Section [5] We do n ot include the CORAVEL data 
(from iLucke fc Mavorlll980t ) since their precision is substan- 
tially worse than ours. 



4 MODELING 

We combined all V 2 and RV measurements in a simultane- 
ous least-squares fit to derive the full orbital solution and 
the physical parameters of BY Dra. We used our own proce- 
dure which minimizes the x 2 function with a least-squares 
Levenberg-Marquardt algorithm. The procedure fits a Ke- 
plerian orbit with corrections to the RVs due to tidal dis- 
tortions of the components and relativistic effects. In or- 
der to mode l the tidal term we use t he Wilson-Devinney 
(WD) code ([Wilson fc Devinnevl Il97ll ) as is explained in 
iKonacki et al. I |2010r i and assume several parameters of BY 
Dra listed in Table [2] Note that both the relativistic and 
tidal effects are much smaller than the RV scatter (see Fig.[T|) 
but we decided to include them in the RV model anyway to 
maint ain a consistent trea tment of our iodine cell based RVs 
as in lKonacki et al.l (|2010T ). Apparent stellar diameters were 
assumed to agree with the estimates of the radii from Sec- 
tion 6.2, since the components are too small and act like 
point sources. 

For a combined V 2 +RV solution our software evaluates 
the the period P, standard Keplerian elements: major semi- 
axis a (of B relatively to A - apparent astrometric in mas), 
inclination i, eccentricity e, longitude of pericenter w, lon- 
gitude of ascending node f2, time of periastron passage T p ; 
velocity amplitudes K\ and K2, systemic velocity vo, flux 
ratios in the observing bands th and tk, and a set of shifts 
in radial velocities between the two components as well as 
between the data sets from each telescope/spectrograph. On 
this basis the software calculates such absolute physical pa- 
rameters like the absolute major semi-axis ai,2 (relatively to 
the baricentre - in AU), absolute components' masses M\ 
and M2, magnitude differences AH and AK, and parallax n. 
The uncertainty of every parameter is a combination of for- 
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Figure 1. Radial velocity variations of the primary (solid lines) 
and secondary (dashed lines) stars of BY Dra (HD234677) as a 
function of the orbital phase. The top panel (a) shows the RV 
variations due to the tidal distortion of the stars and the bottom 
panel (b) due to the combined gravitational redshift and trans- 
verse Doppler effects which together are the dominant term of the 
relativistic correction. 



Table 2. Assumed parameters of BY Dra 



Parameter 


Primary 


Secondary 


Effective temperature, T (K) 


4000 


4000 


Potential, Q 


24.0 


24.5 


Synchronization factor, F 


1.95 


1.95 


Gravity darkening exponent, g 


0.3 


0.3 


Albedo, A 


0.5 


0.5 


Apparent diameter, 8 (mas) 


0.6 


0.5 


Metallicity 




0.0 



mal best-fit least-squares errors a nd systematic errors as is 
explained in lKonacki et al] l|2010h . For the systematic errors 
we assumed the following estimates for additional uncertain- 
ties related to the V 2 data reduction (1) 0.01 percent in the 
baseline vector coordinates, (2) 0.5 percent in A and (3) 10 
percent in the calibrator and binary components diameters. 
For the RVs we assumed (4) 10 percent in all the parameters 
from Table 2 except for the temperatures for which we as- 
sumed an uncertainty of 2 percent and for the metallicities 
we assumed an uncertainty of 0.05 dex. 



5 RESULTS 

The results of our modeling are collected in Tables [3J and 
[4] Figure [2] shows our RVs together with the best fitting 
orbital solution and the corresponding residuals and their 
histograms. Figure [3] shows the same for the PTI V 2 mea- 
surements. The resulting astrometric orbit of component B 
relative to A is shown in the panel (d). In Table |3j we show 
the orbital parameters for BY Dra, the velocity offsets and 
other parameters related to the quality of the fit. The abso- 
lute physical parameters are listed in Table [4] 

As one can see, we were able to reach ~0.2 % of preci- 
sion in velocity amplitudes, despite such obstacles like the 
presence of spots or some rotational broadening of spectral 
lines. This level of quality has direct implication for the pre- 
cision of mass ratio q (0.28 %), Msin 3 i (0.44 and 0.46 % 
for the primary and secondary respectively) or major semi- 
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Figure 2. Observed and modeled radial velocities of BY Dra 
as a function of the orbital phase (a), and their best- fit resid- 
uals as a function of the orbital phase (b) and time (c). 
The histograms of the residuals for the primary and secondary 
(d). The Keck I/HIRES measurements are denoted with cir- 
cles, Shane/CAT/HamSpec with triangles and TNG/SARG with 
stars. Colour version of the Figure is available in the on-line ver- 
sion of the article. 



axis (~0.2 % both for the apparent and absolute values). 
The level of precision in a also proves that the quality of the 
astrometric solution is very high. The 299 visibility mea- 
surements used provide good orbital phase coverage. The 
apparent and physical values of major semi-axis allow us 
to determine the parallax, thus the distance to the system, 
with a precision also close to 0.2 %. Our value of the paral- 
lax - 60.43(12) mas - is in a relatively good agreement but 
almost 6 times more precise than 61. 15(68) mas from th e 
new reduction of the Hipparcos data (jvan LeeuwerJl2007h . 
We were also able to precisely derive the magnitude differ- 
ence in the K band (282 V 2 measurements) but the accuracy 
for the H band is much lower due to a lower number of V 2 
measurements in H (only 17). 

Our final error in the absolute masses of the BY Dra 
components is however much higher - 3.3 % for both the 
primary and secondary. This is mainly due to the inclina- 
tion of the orbit of 154.4 deg. For such configurations, far 
from edge-on, a small error in the angle propagates to a 
large error in the masses. Still it i s a considerably more ac- 
curate measurement compared to iBoden fc Land (|200ll ) of 
respectively 23% and 25% for the primary and secondary. 
This is possible thanks to our superior RV data set (rms of 
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Figure 3. Visibility measurements of BY Dra as a function of 
time (a), and their best- fit residuals as a function of time (b) 
and histogram (c). The 299 measurements used to determine the 
best-fit orbital solution are denoted with black filled circles. The 
corresponding orbital coverage and the relative orbit is shown in 
the panel (d). Colour version of the Figure is available in the 
on-line version of the article. 



~0.15 km s _1 vs 2.3 km s" 1 ) 
astrometric V 2 data. 



6 DISCUSSION 



and a longer time span of the 



6.1 Age and metallicity 

As it is pointed out bv I Torres. Andersen fc Gimened <|2010h . 
the mass uncertainty should be below 3 % to be useful to 
perform reliable tests of the stellar evolution models. Our 
precision is close to that but to go below 3 % we would 
require a higher number of precise RVs or more V 2 mea- 
surements. Nevertheless, with our measurements we still can 
place some constraints on the evolutionary properties of BY 
Dra. We focused on the age estimation to confirm or exclude 
the pre-main-sequence nature of the system. 

We comp ared our results with the Yonse i-Yale 
isochrones (Y 2 ; lYi et al]|200ll ; IPemarque et aT1l2004h . We 
used our estimations of the magnitude differences, paral- 
lax and the appa rent H and K magnitude from 2MASS 
l|Cutri et al.l 120031 ) to derive the absolute H and K mag- 
nitudes of each compon ent separat e ly. Us ing the transfor- 
mation equations from ICarpenterl (|200ll . with updatesjf] 
we transformed them to the ESQ ph otometric system 
l|van der Bliek. Manfroid fc Bouchedll996t ) , in which the Y 2 
isochrones are available. In Figure U we show our measure- 
ments in the mass/A'-band (left) and mass/i/-band (right) 
absolute magnitude diagrams. For comparison we plot the 
isochrones for ages of 60 Myr (dashed) and 1 Gyr (solid 
lines) for three values of the metallicity: Z = 0.02 (red), 
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Figure 4. Comparison of our results with the Yonsei-Yale isochrones in the mass/if- band (left) and mass///- band (right) absolute 
magnitude (in ESO system). The isochrones for 1 Gyr are depicted with solid and for 60 Myr with dashed lines. The isochrones for 
Z = 0.02 are depicted with red, for 0.04 with green and for 0.06 with blue lines. Colour version of the Figure is available in the on-line 
version of the article. 



Table 3. Best-fitting orbital solution and its parameters for BY 
Dra. 



Table 4. Physical parameters of BY Dra 



Parameter 



Value(±) 



Orbital solution 
Apparent major semi-axis, a (mas) 
Period, P (d) 

Time of periastron, T p (TDB-2450000.5) 
Eccentricity, e 

Longitude of periastron, uj (deg) 
Longitude of ascending node, f2 (deg) 
Inclination, i (deg) 

Magnitude difference in K band, AK (mag) 
Magnitude difference in H band, AH (mag) 
Velocity amplitude, primary, K\ (km s _1 ) 
Velocity amplitude, secondary, K2 (km s —1 ) 
Mass ratio, q 

Gamma velocity, vq (km s _1 ) 



4.4472(91) 
5.9751130(46) 
3999.2144(21) 
0.30014(62) 
230.33(17) 
152.30(10) 
154.41(29) 
0.530(11) 
0.60(23) 
28.394(60) 
32.284(61) 
0.8795(25) 
-25.484(46) 



Velocity offsets (all in km s 1 ) 
Secondary vs primary 
SARG vs HIRES, primary 
SARG vs HIRES, secondary 
HamSpec vs HIRES, primary 
HamSpec vs HIRES, secondary 

Least-squares fit parameters 
Number of RV measurements, total 
Number of RV measurements, HIRES 
Number of RV measurements, SARG 
Number of RV measurements, HamSpec 
Number of V 2 measurements 
Combined RV rms, prim. /sec. (km s -1 ) 
Visibilities V 2 rms 
RV x 2 , primary/secondary 
Visibilities V 2 \ 2 
Degrees of freedom, DOF 
Total reduced x 2 , X 2 I DOF 



-0.088(67) 
-0.216(104) 
-0.343(105) 
0.076(83) 

-0.067(85) 



48 
30 
8 
10 
299 

0.169/0.157 
0.0312 
29.17/24.65 
395.3 
330 
1.361 



Parameter 



Primary 



Secondary 



Major semi-axis, a 
Major semi-axis, a 
M sin 3 i (M Q ) 
Mass, M (M Q ) 
M K ,2MASS (mag) 
M H ,2MASS (mag) 
Parallax, k (mas) 
Distance, d (pc) 



(10~ 2 AU) 
(R©) 



3.4437(73) 
7.400(16) 

0.06387(28) 
0.792(26) 
4.269(21) 
4.420(86) 



3.9155(74) 
8.414(16) 

0.05618(26) 
0.697(23) 
4.799(22) 

5.020(149) 



60.43(12) 
16.548(35) 



0.04 (green) and 0.06 (blue). One can see that the prop- 
erties of the secondary component are not reproduced by 
the 60 Myr isochrones, which means that it has already set- 
tled down on the main sequence. The formally best match 
is found for t — 1 Gyr and Z = 0.04. No match was found 
for ages below 60 Myr for any Z value nor for any age value 
for Z < 0.02. For t > 5 Gyr only isochrones with metallic- 
ities higher than solar reproduce the data points. We thus 
conclude that BY Dra is probably between 0.2 and 5 Gyr 
old and is more metal-rich than the Sun. The most probable 
values of age and Z are 1-2 Gyr and 0.04 respectively These 
facts make the pre-main-sequence scenario less probable. 

One should notice that the error bars in the masses 
are enlarged mainly by the uncertainty in the inclination. 
Any change in i would shift both components in the same 
direction - towards higher or lower masses which would def- 
initely not improve the fit. We also have a large uncertainty 
in Mjj, induced by the error in AH, which is so large due to 
a small number of V 2 measurements in this band. Reduction 
of this uncertainty would allow for putting even more strin- 
gent constrains on the nature of the system. At the same 
time, AK is very well constrained and shows that the mass 
ratio q ~ 0.88 is not inconsistent with the observed flux ra- 
tio, at least in the K band. Using the Y 2 set of isochrones 
we can estimate that the expected theoretical magnitude 
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Figure 5. Radii of the components as predicted from our es- 
timations of BY Dra's parameters and the Y 2 best-matching 
isochrone. Horizontal solid lines represent our estimations of the 
masses and the dashed lines the 1-sigma ranges of their respective 
uncertainties. 



difference in V for the stars having 0.792 and 0.697 M5 



1.15±0.1 mae predicted bvlLucke & Mavoi 


J 19801). However 


given even ~0.2 mag variation from spots ( 


ChuEainovll 19661; 



such a difference in V is possible for BY Dra even if it is a 
main-sequence system. 

To put additional constrains on the system's age, 
we further calculated the galactic space velocities 
U, V, W0 relatively to the local standard of rest (LSR; 
I Johnson fc Soderbloml [l987h . We applied our values of ra- 
dial systemic velocity and distance estimation together with 
proper motion of /j, a = 185.92 mas yr -1 and us = —324.8 1 
mas yr" 1 from the PPMX catalogue (|Roser et all l2008h . 
Values of U = 28.2 ± 0.1, V = -13.16 ± 0.06 and W = 
—21.75 ± 0.10 km s _1 put BY Dra out side of any know n 
young moving group or group candidate (|Zhao et al.l f2009). 
and at th e transition area be t ween the thin and thick galac- 
tic disk jBensbv et alj|2003l : [Nordstrom etldl l2004h . This 
supports the possibility of BY Dra being not a PMS sys- 
tem. 



6.2 Spin-orbit (mis)alignment 

Using the masses and isochrones, we can estimate the radii 
of each component of BY Dra. In Figure [5] we plot the Y 2 
isochrone for 1 Gyr and Z — 0.4 in a radius/mass plane. As 
solid horizontal lines we plot the masses together with their 
uncertainties (dashed lines). Other probable isochrones are 
very close to the chosen one and do not change the results 
of our analysis significantly. 

Our results predict values of Ri = 0.695 ± 0.025 and 



R2 ~ 0.61±0.02 R0. We can use this together with the rota- 
tional velocities and rotational period values from the litera- 
ture to estimate the orbital inclination angles. For this pur- 
pose we use v ro t,i sinvot. i = 8.05±0.33 andiy t ,2 sini rotj 2 = 
7.42 ± 1.06 km s" 1 fromlLucke fc Mayor] (Il98fj|j and P ro t = 
3.8285 d from iPettersen. Olah fc Sandmannl l| 19921 ). This 
implies R\ sin i ro t,i = 0.61 ± 0.03 and R2 sini roi ,2 = 0.56 ± 
0.08. Most of the authors attribute spots to the primary 
component and refer the ~3.8 d period to its rotation. If so, 
from the values above we can estimate the rotational incli- 
nation of the primary to be i ro t,i = 6lly 2 or 119^12 deg. 
If we assume the secondary to rotate with the given period, 
we end up with i ro t,2 = 67 or 113 deg and its uncertainties 
ranging from 50 to 130 deg. This means that i ro t,2 = 90 deg 
is also possible . Our results are in a good ag reement with 
values given bv lGIebocki fc Stawikowskil (|l995 ) who derived 



irot,i = 60_ 9 and i ro t,2 = 85_ 15 deg. None of the values 
however agrees with the orbital inclination i oro = 154.4 deg. 
This indicates the spin-orbit misalignment in the BY Dra 
system. Even if we consider that the theoretical radii are 
underestimated by about 15%, a well known issue for late- 
type stars, we are not able to reproduce the observed i OT b- 

In the case of a spin-orbit alignment, and the literature 
values of v ro t sin i ro t, the radii would have to be -Ri ~ 1.41 
and R2 — 1.30 Rq. This would occur if the system was 
~3-4 Myr old, depending on the metallicity. In such a case 
both stars should be much brighter in the infrared than it 
is observed. The predicted if-magnitudes difference for the 
two stars would be ~0.2 and not 0.53 mag which is observed. 
One also would expect the system to be a member of a young 
cluster, containing leftovers of the primodial gas, but this is 
not observed as well. 

However, the spin-orbit alignment sho uld be ob served 
for such a close pair of ~1 Gyr old stars (|Hutlll98lh . The 
source of the discrepancy could be for example an overes - 
timated rotational velocity. iGIebocki fc Stawikowskil (|l995l ) 
in their analysis a dopte d a value of 3.6 km s _1 , given by 
IStrassmeier et all (|l993l ri For this value of v ro t sin i rot we 
get i ro t,i — 157 deg, which is very close to the observed or- 
bital inclination. Hence the spin-orbi t alignment may in fact 
be present in the BY Dra system if iLucke fc Mayor! (Il980l ) 
have overestimated their rotational velocity measurements 
by adopting a too small macroturbulence velocity. 

Finally, let us note that a spin-orbit misalignment could 
manifest itself throu gh its impact on the apsidal precession 
rate (for a review see lMazehll200sl) . Unfortunately since BY 
Dra is not an eclipsing system and our RVs and V 2 s are 
not sufficiently accurate, a measurement of the apsidal mo- 
tion cannot be carried out. We attempted to fit for uj but 
obtained statistically insignificant value. 



6.3 Rotation pseudo-synchronization 

In the case of eccentric orbits one can find a rotational pe- 
riod for which an equilibrium is achieved. This equilibrium, 
called the pseudo-synchronization, occurs when the ratio of 
the orbital to the rotational period is: 



4 Positive values of U, V and W indicate velocities toward the 
Galactic center, direction of rotation and north pole respectively 



5 In fact, in t he c urrent version of t heir catalogue, 
Strassmeicr ct al. cite v ro t values from ILucke &: Mayor] 

jl980h . 
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1 + 7.5e 2 



5.625e 4 + 0.3125e 6 



Prat, P s (l-e 2 )3/2(l + 3e 2 + 0.375e 4 ) (5) 

(|Hutl Il98ll ; iMazehl 120081 ). For the observed eccentricity 
of BY Dra we get P or b/Prot, P s = 1.559(3) or P r ot, ps = 
3.833(8) d, which is in a good agreement with the observed 
Prot = 3.8285 d, and Porb/Prot = 1.561. One may thus 
conclude that BY Dra is in a rotational equilibrium. The 
predicted time scale of the pseudo-synchronization is in the 
case of BY Dra s imilar but slightly shorter than for the spin- 
orbit alignment (|Hutlll981|) ■ Using the approximate formula 
for the synchroniza tion time scale of late-type stars given by 
iDevor et all (^POSf ), we get the value of the order of 10 Myr, 
so still smaller than the age of the system. The above equa- 
tion was however derived for binaries with no additional 
companions (see below) and the tidal evolution of BY Dra 
AB might be different if the gravitational influence of the 
third body is taken into account. 



6.4 Eccentricity and multiplicity of BY Dra 

The eccentricity of BY Dra AB (e = 0.3) appears to be 
unusually high. According to lZahn fc Bouchetl |l989) a cir- 
cularization of the orbit of a late type system such as BY Dra 
should occur during the pre-main-sequence phase. This was 
on e of the argument s for t he PMS nature of BY Dra. Based 
on IZahn fc Bouchetl £l989) we can estimate that in the case 
of BY Dra the eccentricity should drop to a few percent over 
~ 10 5 years. 

However, BY Dra is a hierarchical triple system, with 
a distant common proper motion companion. The projected 
separation of 16.7 arc-sec and our distance determination 
indicate a projected physica l separ ation of 277 AU. As esti- 
mated by IZuckerman et all ([1997) its mass is about 0.13 
Mq and assuming a circular orbit, it corresponds to an 
orbital period of about 2050 years. It is conceivable that 
the observed eccentricity of the BY Dra AB pair could be 
explained by the presence of the companion through the 
Mazeh-Shaham mechanism which results in a cyclic eccen- 
tricity variation, known as the Kozai cycles (|Kozai| 1962; 
Mazeh fc Shahamll 19791 ; iFabrvckv fc Tremainell2007l ; IMazehl 



K x = 157 m/s 



K x = 78 m/s 



2008f ). 



Let us denote all the orbital and physical parameters 
of an unknown perturber by the index X. In order to put 
some constrains on the properties of the perturber which 
would induce sufficiently strong Kozai cyc les, we followed 
the analysis of IFabrvckv fc Tremainel (j2007). The two main 
conditions which have to be met in order to produce the ob- 
served eccentricity are: (1) a sufficiently large relative incli- 
nation i re i of the binary (inner) and the perturber's (outer) 
orbit; (2) the Kozai cycles time-scale r must be shorter 
than the period of the inner orbit's pericenter precession. 
For the BY Dra AB pair the relativistic precession is the 
dominant one, being at least an ord er of magnitude faster 
than any other (tidal or rotational; IFabrvckv fc Tremainel 
l2007h . The precession period is Ugr — 25000 yr, which cor- 
respon ds to uigr = 8.0 x 10~ 1 2 rad s _1 . Using the formal- 
ism of IFabrvckv fc Tremainel ([20071 ) we can estimate that 
the observed eccentricity of BY Dra can be induced when 
TWGfl|e=o ^ 2.796 (in SI units), where the term ojgr is com- 
puted for e = 0. From this we can derive the following con- 
ditions for the parameters of the perturbing body: 



1000 2000 

p x [d] 
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Figure 6. Conditions necessary to reproduce the observed eccen- 
tricity of BY Dra through the Mazeh-Shaham mechanism by bod- 
ies which would produce the RV signal with the semi-amplitude 
of 157 and 78 m s — 1 (the lower rms of the RV orbital solution and 
half of its value, left and right column respectively), and having 
orbital inclinations of 90 (top row) and 52.4° (bottom row). The 
shaded areas refer to the periods and eccentricities which would 
not allow to induce sufficiently strong Kozai cycles. The lower 
limit of the period is 35 d (the semimajor axi s of 0.24 AU) which 
refers to the shortest stable circular orbit l lHolman fc Wiegertl 
1999). The area above the solid line corresponds to the eccentric 
orbits whose periastron distance is within the instability zone (i.e. 
is shorter than 0.24 AU). The upper limit in the eccentricities 
shown in the figures is 0.9. 



2 Mi + M 2 + Mx 
Of - 



Mx 



e 2 x ) 3/2 < 938.21 



u x 
M x 



ex 



3/2 



< 938.21, 



(6) 



(7) 



where the orbital period is given in years, major semi-axis 
in AU and all masses in Mq. The condition rcjG_R| e= o ^ 
2.796 also allows us to deduce that the relative inclina- 
tion of the two orbits must be larger than 78° (or smaller 
than 102°). This value c an be confirmed by the results of 
iFord. Kozinskv fc Rasiol l|2000h . which for the mass ratio of 
BY Dra AB predict i re i > 75°. The relatively narrow range 
of i re i allows us to put some usefull constrains (which in- 
clude possible values of the longitudes of ascending nodes) 
on the absolute value of the perturber's orbital inclination: 
ix £ [52.4°, 127.6°]. 

In Figure[5]we present results of our analysis. We use re- 
lation[6]to check whether a body of a given orbital properties 
(i.e. period and eccentricity) can produce sufficient Kozai cy- 
cles. For the mass Mx we have taken mass of a body that 
in an orbit of a given Px, ex, and ix = 52.4 or 90° would 
produce an RV modulation of the inner pair at the level of 
the smaller rms (157 m s _1 ) or half of that. The four pan- 
els show the Px / ex parameter space for the two values of 
inclinations and RV semi-amplitudes. The shaded areas cor- 
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Figure 7. Limits for Mx sin(ix) for a putative fourth body es- 
timated from our RV measurements, for several values of eccen- 
tricities. The value of the rms = 157 m s _1 was assumed to 
be the RV semi-amplitude. The lower limit of the period is 35 
d (the semimajor ax is of 0.24 AU) which refe rs to the shortest 
stable circular orbit jHolman fc Wieeer j[l99g|) . For eccentric or- 
bits the limits are terminated at the shortest periods having the 
distance of pcriastron larger than 0.24 AU (around 60, 139 and 
1107 days). Using ix = 52.4° increases the detection limits by 
about 26% (l./sin(52.4°) ~ 1.26). Colour version of the Figure 
is available in the on-line version of the manuscript. 

respond to the values of Px and ex for which the observed 
eccentricity of BY Dra AB would not be induced. The solid 
line shows the short-period stability border, calculated in 
such way that for a given eccentricity the orbit has its peri- 
astron at ~0.24 AU which r efers to the smallest stab le cir- 
cular orbit (with P ~ 35 d; iHolman fc Wiegerdll999h . The 
long-period cut-off at 2900 days is mainly for the clarity, but 
it is close to double time span of our observations (1513 d) 
which means that RV modulations with periods around 3000 
d could in principle be detected. The corresponding fourth 
body detection limits in terms of Mx sin(ix) as a function 
of its orbital period and eccentricity are presented in the 
Figure 

From those two figures one can deduce that if the ob- 
served eccentricity is an effect of the Mazeh-Shaham mech- 
anism, the perturber should be in an orbit of up to single 
years (major semi axes from ~ 0.2 to ~ 2 AU) and have 
its mass in the planetary regime. This is consistent with 
the fact that the interferometric V 2 measurements are fully- 
consistent with a two-disc model, so no additional light is 
detected. It is not however excluded that the eccentricity 
of the perturber is very large which would allow more mas- 
sive bodies in long period orbits to induce the Kozai cycles 
and remai n undetectable by t he RV s. The component C dis- 
covered bv lZuckerman et al.l |l997) could be the perturbing 
body but with the relation [7] it seems that it is not very 
likely. For the estimated mass M3 = 0.13 Mq the orbital 
eccentricity ez would have to be larger than 0.98, assuming 
a fortunate but improbable case that the star is currently 
seen exactly at the apocenter, and Q.3 = 0°, thus 03 [AU] 
= 277/(l + e3). For less fortunate cases the value of e$ would 
have to be even larger. 

The other possibility is the presence of a putative fourth 



body reported in the Hipparcos catalog. iBoden fc Land 
I2OOII ) inspected the available RV data in order to find a 114 
day period predicted by the Hipparcos catalog and with a 
high level of confidence they excluded the existence of such a 
period in the spectroscopic data. With the relation [6] we can 
show that if a body in a Ph = 114.02 d circular orbit exists, 
it would have to have Mh > 0.15 Mjr/p. From the Figure [7] 
we can put an upper mass limit of 5.44 Mjup, taking into 
account the reported inclination of 113.21°, which itself is 
within the allowed limits. The reported major semi-axis of 
the Hipparcos photocentric orbital solution is ai2,H = 0.0515 
AU (at the distance to BY Dra). Assuming the maximum 
mass ratio Mh/(Mi+M 2 ) = 0.0035, the 4-th body barycen- 
tric major semi-axis would be 16.07 AU (after correcting for 
the inclination), but the Kepler's 3'rd law predicts the ma- 
jor semi-axis of 0.53 AU for the given period and masses. 
Such a body would produce the RV signal much stronger 
than 157 m s _1 . We can thus conclude that the Hipparcos 
solution is unrealistic. 



7 SUMMARY 

We present the most precise orbital and physical param- 
eters of an important astrophysical object - the low-mass 
SB2 BY Draconis, a prototype of an entire class of vari- 
able stars. We reach a level of precision which allows us to 
put important constrains on the nature of this object. We 
conclude that this is not a pre-main-sequence system, de- 
spite its high orbital eccentricity and a possible spin-orbit 
misalignment. However, the gravitational influence of a 4-th, 
yet undetected body in the system may explain the observed 
value of e and the spin-orbit alignment may be inferred from 
the available data if a smaller than claimed in the literature 
value of the rotational velocity is used. The observed rota- 
tional period and the eccentricity suggest that the BY Dra 
AB system is in the rotational equilibrium. However, if the 
observed eccentricity is indeed due to the presence of the 4- 
th body, the putatiuve companion may have its mass in the 
planetary regime. The whole dynamical and tidal picture of 
this system is more complicated than we previoulsy thought 
and deserving perhaps a dedicated theoretical, observational 
and numerical analysis. 
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cellent condition and operating PTI. 
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